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Abstract. The evolution of a zero metallicity 9 Mq star is computed, analyzed and compared with that of a solar 
metallicity star of identical ZAMS mass. Our computations range from the main sequence until the formation of 
a massive oxygen-neon white dwarf. Special attention has been payed to carbon burning in conditions of partial 
degeneracy as well as to the subsequent thermally pulsing Super- AGB phase. The latter develops in a fashion 
very similar to that of a solar metallicity 9 Mq star, as a consequence of the significant enrichment in metals 
of the stellar envelope that ensues due to the so-called third dredge-up episode. The abundances in mass of 
the main isotopes in the final ONe core resulting from the evolution are X(^®0)~ 0.59, X(^''Ne)« 0.28 and 
X(2*Mg)« 0.05. This core is surrounded by a 0.05 Mq buffer mainly composed of carbon and oxygen, and on top 
of it a He envelope of mass ~ 10~* Mq. 
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1. Introduction 

. The primordial universe was mainly composed of hydrogen 
; {X « 0.77), helium {Y « 0.23), and some ^H, ^He and ^Li, 
. in a very small amount — the total abundance did not rep- 
'■ resent more than 10~^° of the total mass of the universe. 
. This set of abundances, derived from the standard ho- 
' mogeneous Big Bang nucleosynthesis (Alpher & Herman, 
1950; Olive, 1999), characterizes the matter out of which 
the first structures of the Universe formed. According to 
the accepted hierarchical scenario, these large structures 
had a mass between 10^ and 10* Mq, and very soon after 
their formation became the nest in which primordial stars 
formed. 

To this regard it is important to mention at this point 
that there is no general agreement on the initial mass func- 
tion (IMF) of these primordial stars. Instead, this par- 
ticular subject is still today a matter of strong debate. 
Moreover, the original composition of these stars was def- 
initely determined by the fragmentation and condensation 
processes that led to their formation and, therefore, it is 
closely related to the adopted IMF. The first attempts 
on this direction pointed to the dominance of very mas- 
sive stars. The reason was the original lack of metallic 
molecules that allowed the cooling necessary to form low 
and intermediate mass objects. Therefore, a big deal of 
the earliest theoretical work on Z = stars was devoted 
to the study of very massive objects — see, for instance, 



Ezer & Cameron (1971) and Carr et al. (1983). However, 
as it was claimed later, the mere presence of molecular 
hydrogen could allow the cooling necessary for low and 
intermediate mass objects to be formed (Carlberg, 1981; 
Palla et al., 1983; Yoshii & Saio, 1986). This result moti- 
vated the work of D'Antona (1982), Guenter & Demarque 
(1983) and Fujimoto et al. (1984) on intermediate mass 
primordial stars. The latest two-dimensional simulations 
(Nakamura & Umemura, 2001), have found that a bimodal 
IMF peaked at ~ 2 Mq and at a mass between 10 and 
100 Mq could be appropriate. Consequently, the most re- 
cent publications keep studying both massive (Heger et 
al., 2001; Heger & Woosley, 2000) and low- intermediate 
mass objects (Marigo et al., 2001; ChiefR et al., 2001). For 
instance, Heger et al. (2000) considered the evolution of 
Z — stars of masses between 15 and 250 Mq, from their 
main sequence until core collapse. They concluded that 
massive objects in general were important in the produc- 
tion of primordial ^"^N and that those stars within the 
highest mass range were the place were primordial Fe and 
Ni could be synthetized — see also Goriely & Siess (2002) 
and Abia et al. (2001). Using the results of Limongi et al. 
(2000) and Chieffi et al. (2001) for Population HI stars, 
Abia et al. (2001) analyzed the impact of primordial stars 
in polluting the intergalactic medium (IGM), and showed 
the relevance of an IMF peaking in the range of masses 
between 4 and 8Mq in order to explain the large [C/Fe] 
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and [N/Fe] ratios observed in many extremely metal-poor 
stars. 

In a pioneering work on low and intermediate- mass 
stars of very low metallicity, Fujimoto et al. (1984), stud- 
ied by means of semianalytical techniques the generation 
of helium-shell flashes during the AGB stage, as a func- 
tion of mass and metallicity. They found that stars whose 
ZAMS mass was lower than 4 Mq did experience helium 
shell flashes and were important contributors of ""^^C and 
neutron-rich isotopes of light elements, whereas stars with 
larger masses were more likely to become supernovae and, 
hence, to contribute substantially to the production of Fe 
and Ni. Marigo et al. (2001) and Marigo (2002) also stud- 
ied Z = stars in the mass range 0.7 ^ M/Mq ^ 100, 
and followed their evolution from the main sequence until 
the AGB phase for objects in the low- and intermediate- 
mass range. In particular, they considered the pollution 
of stellar atmospheres and, ultimatelly, of the interstel- 
lar medium, as a consequence of the dredge-up processes 
ensuing during the AGB phase. For the high-mass range, 
namely stars with masses Mzams <^ 9-^-^0; Marigo et al. 
(2002) followed the evolution of their model stars until the 
carbon burning phase. It is, however, important to stress 
at this point that in the high-mass range carbon is burnt 
in non-degenerate conditions. Chieffi et al. (2001) focused 
on Z = 0, 4 ^ M/Mq ^ 8 stars and also followed the 
thermally pulsing phase and the dredge-up during the 
AGB phase. These authors concluded that stars in this 
mass range were main contributors of nitrogen and oxy- 
gen to the IGM. Finally, Siess et al. (2002) also followed 
the evolution and nucleosynthesis of primordial stars of 
masses ranging between 0.8 and 20 Mq. Their results, in 
agreement with those of Chieffi et al. (2001), pointed to a 
standard AGB behaviour for stars with masses larger than 
about 5 Mq. These authors have also discussed the effects 
of taking into account the effects of a moderate amount 
of overshooting. 

Finally, it is important to point out that one might 
argue the scarcity of observational support against the 
study of Z = objects. In fact, the detection of primordial 
composition objects should take us back to [Fe/H]?s —8.3 
epochs, whereas, up to now, the oldest stars detected are 
a red giant at [Fe/H]« -4.6 (Bessel & Norris, 1984), 
and a low mass star, HE0107-5240, with [Fc/H]w -5.3 
(Christlicb ct al., 2002). Nevertheless, as correctly pointed 
out by Chieffi et al. (2001), the lack of observations should 
not let us discard the existence of Z = intermediate mass 
stars. Therefore, in this paper we study the evolution of 
intermediate mass primordial stars. 

In order to keep consistency with our previous re- 
sults, we have used for the calculations reported in this 
paper the same evolutionary code as that described in 
Ritossa, Garcia-Berro & Iben (1995). The interested 
reader can find there a thorough description of the most 
relevant physical inputs (neutrino energy loss rates, reac- 
tion rates. . . ). In particular, the treatment of convective 
boundaries turns out to be important for the calculations 
reported here. We have not adopted any overshoot and 



convective regions were computed according to the stan- 
dard procedure given by the mixing length theory in order 
to adequately compare with the c;alculations reported in 
Ritossa et al. (1995). The only major difference between 
the physics used there and those used in the present calcu- 
lations is that, in order to properly compute the evolution 
of Z = model star sequences, we have implemented the 
corresponding opacities at very low metallicities of Iglesias 
& Rogers (1993). 

The plan of this paper is the following: in section 2, we 
study the evolution of a zero metallicity 9 Mq model star 
until the beginning of carbon burning, and we compare it 
with the evolution of the corresponding solar metallicity 
model star and with the results of previous calculations. In 
section 3, the carbon burning phase, the associated ther- 
monuclear flashes, the dredge-up of nuclear processed el- 
ements, and the associated changes in the surface compo- 
sition of the star arc studied. Section 4 is devoted to the 
study of the thermally-pulsing Super-AGB (TP-SAGB) 
phase. In section 5, the evolution through the post- AGB 
track and the final phases of the evolution to become a 
white dwarf are described. Finally, section 6 is devoted to 
summarize our main results and to outline the conclusions 
derived from them. 

2. Hydrogen and helium burning 

We consider a zero metallicity 9 Mq star. As we will see 
below, the characteristics of the evolution during the main 
hydrogen and helium burning phases show important dif- 
ferences with respect to the evolution of analogous cases 
of higher metallicity. In particular, as already noted by 
other authors that have studied the main sequence phase 
of intermediate mass primordial stars — namely Chieffi & 
Tornambe (1984), Tornambe & Chieffi (1986) and Cassisi 
& Castellani (1993) — the first ascent to the red giant 
branch does not take place before central heliTim ignition 
but, instead, it occurs when helium burning is set in a 
shell. 

The main characteristic oi Z = stars, that is, the 
absence of metals and, in particular, of carbon, nitrogen 
and oxygen, determines the evolution from its very early 
stages. This composition does not allow the CNO-cycle 
to ensue at the beginning of the core hydrogen burning 
phase, as it is the case for intermediate mass stars of higher 
metallicity. Instead, the pp-chains must provide the nu- 
clear energy for the star to survive approximately during 
the first 7.4 x 10® yr. However, the pp chains are far from 
being the main nuclear energy suppliers during the bulk of 
core hydrogen burning phase and, as soon as the temper- 
ature threshold for helium burning reactions is reached, 
the production of ^^C, ^"^N, and ^^0 in small amounts (at 
abundance thresholds of about 10~^° in mass), is enough 
for the CNO cycle to take over as the primary energy 
supplier during this phase. It derives from our calcula- 
tions that the entire core hydrogen burning phase lasts 
for about 2.1 x 10^ yr, in good agreement with the result 
of Marigo et al. (2001), who obtained 2.2 x lO'^ yr. 
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Table 1. Characteristics of models at various points in the H-R diagram of Fig. 2 
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Model 




logL 


logTeft 


logi? 


log Pc 


logTc 


Mhc 


Xc(He) 


Mc 


Xc(C) 


A 


0.0000 


3.63 


4.55 


0.25 


4.72 


7.71 


0.00 


0.23 


0.00 


0.00 


B 


2.3545 


3.74 


4.61 


0.17 


5.16 


7.87 


1.90 


0.41 


0.00 


1 X 10"^° 


C 


6.8994 


3.95 


4.53 


0.44 


5.24 


7.93 


1.90 


1.00 


0.00 


3 X 10"^ 


D 


7.0477 


4.03 


4.50 


0.54 


6.63 


8.07 


1.90 


0.99 


0.00 


5 X 10~^ 


E 


7.5843 


4.10 


4.53 


0.55 


6.34 


8.22 


2.03 


0.48 


0.52 


0.45 
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8.0000 


4.15 


4.47 


0.65 


7.05 


8.40 


2.10 


0.00 


0.98 


0.36 
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8.0520 


4.08 


4.14 


1.29 


8.80 
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2.15 
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1.05 


0.36 
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Fig. 1. Evolution in the logpc — logTc plane of the zero 
metallicity 9 Mq star (solid line) and oiz,Z = 0.02 star of 
the same mass (dashed line), up to the point where carbon 
is ignited. 



The most important parameters relevant for a cor- 
rect understanding of the structure and composition of 
our models during the hydrogen burning phase are shown 
in Table 1. For the sake of completeness we also show 
in Tabic 1 the same parameters for the rest of the evo- 
lutionary phases before carbon is ignited. The positions 
of each labelled model are also shown in Fig. 1 and Fig. 
2, wliic;li display, respectively, the evolution of the central 
temperature as a function of the central density for both 
our zero metallicity 9 Mq star and that of a Z = 0.02 
star of the same mass, and the corresponding evolution in 
the Hertzsprung-Russell diagram for each of these stars. 
Models labelled as A, B and C correspond, respectively, 
to the beginning of the core hydrogen burning (hereafter 
CHB), to the time at which the energy production through 
the CNO-cycle starts to take over that of the pp-chains 
and, finally, to the end of CHB. Models labelled as D, 
E and F correspond, respectively, to the beginning, to an 
intermediate stage, and to the end of central helium burn- 
ing (hereafter CHcB). Finally, model G corresponds to the 
time at which carbon is ultimately ignited. 
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Fig. 2. Evolutionary tracks in the Hertzsprung-Russell di- 
agram of the 9 Mq zero metallicity star (top panel) and of 
a star of 9 Mq with Z = 0.02 (bottom panel). The labels 
in the Z = model correspond to the models of Table 1. 
The labels in the Z = 0.02 model correspond to equivalent 
evolutionary stages. 

As it can be seen in Fig. 1, during the phase in which 
the pp-chains are the dominant reactions — from A to 
B — core contraction proceeds accompanied by a smooth 
increase in the central temperature. However, once the 
CNO-cycle is activated (at B) , the larger release of energy 
that consequently ensues temporarily stops the contrac- 
tion and lets the core expand before the usual behaviour of 
core contraction and heating during the bulk of the CHB 
phase is attained. For the sake of comparison, we have 
also represented in Fig. 1 the evolution of the central tem- 
perature as a function of the central density for the case 
Z = 0.02. As it can be seen, important differences arise. 
The loci in the log p — log T diagram where core hydrogen 
burning begins and fully develops are different because 
hydrogen burning through the pp-chains requires higher 
temperatures than through the CNO cycle to achieve the 
same rate of energy production. This fact also explains 
the positions of the associated loops in both curves. Both 
tracks are practically coincident at the point where the 3a 
reaction starts (D), as central He burning is not sensitive 
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Fig. 3. Evolution of the convective regions of our model 
star (bottom panel) and of the radius of the star (top 
panel) as a function of time during the CHB and CHeB 
phases. 
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Fig. 4. Evolution of the central abundances of the 9 Mq 
primary until the beginning of carbon burning. 



to metallicity effects. Nevertheless the subsequent evolu- 
tion takes the Z = star back to the region of higher 
densities. 

The evolution in the Hertzsprung-Russell diagram of 
the zero metallicity star also presents substantial differ- 
ences with respect to that of solar metallicity (see Fig. 2) . 
The part of the core hydrogen burning phase that pro- 
ceeds through the pp-chains appears to occur at slightly 
higher effective temperatures and luminosities than those 
of the Z = 0.02 case. The bulk of the core hydrogen burn- 
ing phase, occuring through the CNO-cycle, leads the star 
to a track whose shape is very similar to that found for 
the Z = 0.02 case. Nonetheless, once this process is com- 
pleted, the star docs not evolve yet to the low Tofr low L 
regions of the Hertzsprung-Russell diagram to reach re- 
gions close to the Hayashi track. Instead, because helium 
burning occurs very shortly after hydrogen is exhausted in 
the core, the star keeps its luminosity almost unchanged, 
and just a small decrease in effective temperature is ob- 
served. The ascent along the red giant branch that can 
be seen in Fig. 2 corresponds to the helium-shell burning 
phase when, as it will be shown in the next section, carbon 
burning has already developed in the core. The former re- 
sults are in agreement with, for instance, those of Marigo 
et al. (2001) or Chieffi et al. (2001). 

The evolution of the surface radius and convective re- 
gions are shown in the top and bottom panels of Fig. 3, 
respectively. The part of the CHB phase during which the 
pp chains dominate can be clearly distinguished from that 
at which the CNO-cycle takes over, as the size of the con- 
vective core is considerably smaller when the pp-chains 
dominate. The pronounced decrease of the size of the cen- 
tral convective region occuring at t « x 10^'* s corresponds 



to a time shortly after the reactions of the CNO cycle set 
in. Moreover, during the phase of CHB dominated by the 
pp-chains {t « 1.45 x 10^^ s) the radius of the star re- 
mains constant, and slightly increases afterwards. Note, 
however, that the surface radius of our model star does 
not increase significantly until CHeB is over. 

The evolution of the central abundances of the main 
isotopes (H, He, C and O) as a function of time up to 
carbon ignition are shown in Fig. 4. The step present in 
the evolution of the central hydrogen and helium abun- 
dances between t m 1.5 X 10^^ s and 2.0 x 10"^^ s is due 
to the competition between two effects. On one hand, hy- 
drogen burning tends to decrease its abundance. On the 
other hand, the very fast increase in the size of central 
convective region that accompanies the activation of the 
reactions of the CNO-cycle allows fresh hydrogen to be 
engulfed and, eventually, keeps its abundance constant. 
The slopes of X(^B.) and X(''He) as a function of time 
are also different during the phase in which the pp-chains 
dominate and during the CNO cycle burning phase. At 
the end of the helium burning phase the He-exhausted 
core has a mass of about 1.1 Mq, and is composed of car- 
bon (X(i2C) ~ 0.34) and oxygen {X(^'^0) :i! 0.65). These 
abundances differ from those obtained by Marigo (2001), 
whose calculations yielded X(i2c)/X(if'0) = 0.658. Also 
the size of the He- exhausted core shown here is smaller 
than that obtained by Marigo (2001), who obtained a 
1.29 Mq He-exhausted core. The differences in the CO 
core masses and abundance ratios are probably a conse- 
quence of the different treatments adopted for convection. 
Marigo (2001) considered a moderate amount of over- 
shooting in their calculations, whereas we are not using 
overshooting at all. 
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Fig. 5. Temporal evolution of the main structural param- 
eters during the bulk of the carbon burning phase. In 

the top panel the luminosity associated to helium burning 
(solid line) and carbon burning (dashed line) are shown. 
The central panel displays the evolution of the stellar sur- 
face radius (solid line) and of the luminosity (dashed line). 
In the bottom panel the evolution of the convective regions 
is represented. See text for details. 

3. The carbon burning phase 

Figure 5 shows the temporal evolution of main structural 
quantities just after the CHcB phase and during the bulk 
of the carbon phase. In particular, in the top panel of 
Fig. 6 the helium and carbon luminosities are shown (solid 
and dashed lines, respectively). The central panel of Fig. 5 
shows the variations of the radius and the luminosity. 
Finally, the bottom panel of Fig. 5 shows the behavior of 
the convective regions during this phase. The vertical thin 
line in all three panels indicates when the carbon burning 
luminosity is 10% of the helium burning luminosity, corre- 
sponding to t = 8.0562 X IQi'' s. Re gions to the left of this 
line correspond to the phase in which the helium burn- 
ing shell is the dominant nuclear energy source, whereas 
regions to the right of this line correspond to phases in 
which carbon burning is significant. 

Once the CHeB phase has been completed, the He- 
burning reactions set in a shell and, simultaneously, 
the star begins an excursion to the red part of the 
Hertzsprung-Russell diagram at slightly decreasing lumi- 
nosities, as shown in Fig. 2 and in the central panel of 
Fig. 5. At this point {t = 8.0571 x lO" s) the surface ra- 
dius of the star has increased from an initial value of about 
20 B.Q up to roughly 320 Rq . This overall expansion is 
accompanied by the formation of a convective envelope, 
whose inner edge is able to penetrate relatively deep in- 
side the star sec the bottom panel of Fig. 5. However, 
the bottom of the convective envelope stops its advance 



at t = 8.0571 X 10^'' s, before it can reach the regions of 
the star that have been previously enriched in products 
of the H nuclear reactions. During this process, when the 
maximum extent inwards of the base of the convective en- 
velope (BCE from now on) is reached (A'/bce = 2.15 Mq), 
the surface abundances (by mass) of carbon, nitrogen 
and oxygen are, respectively, X(^^C) = 8.74 x 10~^^, 
X(14n) = 6.61 X 10-1° and xl^'^O) = 7.84 x lO'^^ 
Therefore, during the early phase of this (so-called) sec- 
ond dredge-up, no significant changes can be observed in 
the surface composition of the star. It must be noted that, 
strictly speaking, we do not find a first dredge up during 
hydrogen shell burning, as it is the case of higher metal- 
licity objects, and consequently it must be stressed that 
calling second dredge-up to this penetration of the con- 
vective envelope is just a matter of keeping consistency 
with the nomenclature used for stars of larger metallicity. 

Carbon burning begins just before the second dredge- 
up, at t = 8.0534 x 10^^ s, while the star is still climbing 
the red giant branch. Note that — as previously men- 
tioned — at time t = 8.0571 x 10^^ s, the surface radius 
stops increasing and the convective envelope temporar- 
ily halts its inner advance. This is due to the increasing 
energy flux released by the vigorous c;arbon burning re- 
actions occuring deep inside the He-exhausted core. As 
a consequence of the ignition of carbon the He-exhausted 
core stops contracting and releasing gravitational energy 
and, consequently, the BCE stops moving inwards. Our 
model star behaves in many senses in a way similar to 
that of solar metallicity stars of slightly higher mass. In 
particular, carbon is burnt off-center in conditions of par- 
tial degeneracy in a series of three consecutive flashes — 
top panel of Fig. 5. The first flash starts at ~ 0.33 Mq 
— bottom panel of Fig. 5 — and the carbon burning lu- 
minosity reaches values as high as Lq — 10^ Lq. 

Fast injections of energy at very high rates and in rela- 
tively small regions of the stellar interior induce high tem- 
perature gradients and, consequently, several inner con- 
vective regions associated to each flash develop — see the 
bottom panel of Fig. 5. Part of the energy released in 
these flashes is transformed into work of expansion that 
changes the physical structure of the layers where helium 
burning is still active. The expansion and cooling of these 
layers lead to an eff'ective decrease in helium luminosity, 
that recovers once the carbon flash is finished, very much 
in the same way solar metallicity stars of the same mass 
do. Nevertheless these dramatic changes affect mainly the 
innermost regions of the star, and they have a rather lim- 
ited effect on the structure of the outer envelope. In fact, 
the surface radius and luminosity remain almost constant, 
similarly to what happens to the Z — 0.02 stars of 10 Mq, 
10.5 Mq and 11 Mq discussed in Ritossa et al. (1995), 
Iben et al. (1997) and Garcfa-Berro et al. (1999), and in 
sharp contrast with the behavior of the 9 Mq star of solar 
metallicity discussed in Garcia-Berro et al. (1997). The 
reason of this behavior is that in the case of the solar 
metallicity star carbon is ignited at Mr ~ 0.46 Mq, much 
closer to the helium-carbon discontinuity than in the case 
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Fig. 6. Relevant structural and dynamical quantities for 
a model at time t = 8.0616 x 10^^ s, just before the second 

advance of the convective envelope, in a region extending 
from below the helium burning shell until the base of the 
convective envelope. Upper panel: velocity (solid line) and 
helium profile (dotted line). Middle panel: time variation 
of the temperature (solid line) and density (dotted line). 
Lower panel: energy release rates (nuclear and gravother- 
mal) along with the time derivative of the internal energy 
and the work of expansion. The nuclear energy release rate 
has been divided by 10 in order to fit into the scale. 



of our Z = model. Hence, for the solar metallicity star 
the carbon flashes are felt all the way to the base of the 
convective envelope much more conspicuously than in the 
Z = model star. 

As a consequence of the second flash, at time t = 
8.06021 X 10^^ s, two important facts occur. First, the 
carbon burning front reaches the center of the star. It is 
worth noting at this point that although the main car- 
bon burning phase (corresponding to the second flash) 
starts indeed as a flash it proceeds as a flame, as dis- 
cussed in depth in Garcia-Berro et al. (1997). Second, 
once the center is reached, a new inward advance of the 
convective envelope takes place (see Fig. 5). During the 
first phase of the dredge-up episode, the base of the con- 
vective envelope reached a mass coordinate of 2.15 Mq, at 
which it stabilized for most of the carbon burning phase. 
After the carbon burning front reaches the center the BCE 
moves down to the regions located deeper in the star. As 
it can be seen in Fig. 5, this penetration of the convec- 
tive envelope occurs in two stages, the first one reaching 
Mr = 1.74 Mq, where it stops temporarily, and the second 
one reaching = 1.20 Mq. In fact the BCE gets deep 
enough to reach regions of the star where hydrogen and 
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Fig. 7. Same as Fig. 8, but for time t = 8.0629 x 10^'^ s, 
during the second inner advance of the envelope. 



helium burning have efficiently proceeded and, therefore, 
the elements of the CNO-cycle and the a-products are 
dredged-up to the stellar surface. The final surface abun- 
dances of carbon, oxygen and nitrogen after this additional 
dredge-up episode are, respectively, X(^^C)= 2.04 x 10"'', 
X("N)= 2.41 X 10-*^ and X(^^0)= 3.47 x 10-'^. 

As already discussed, the first inwards advance of the 
BCE took place during the helium burning shell phase. 
To reach a better understanding on how the dredge-up 
process occurs, wc have represented in Figs. 6 and 7 the 
luminosity and helium abundance profiles (upper panel), 
the relative variation of the temperature and density with 
time (middle panel), and the different rates of energy 
generation and absorption (bottom panel). Fig. 6 corre- 
sponds to a time t = 8.0616 x 10^^ s, just after the car- 
bon burning front has reached the center of the star. As 
it can be seen, the generation of nuclear energy only oc- 
curs at the helium burning shell (HeBS), which is located 
at Mr ~ 1.1926 Mq. This energy is partially absorbed 
in the region above the HeBS (of thickness 0.2 Mq) and 
produces a slight expansion of these layers. Between this 
zone and the BCE, the structure variables T and p and 
the luminosity remain approximately constant. 

Fig. 7 shows the same quantities as in Fig. 8, but at a 
time t — 8.0629 x lO"'^ s, during the second inner advance 
of the envelope. The situation is now very different since 
at this moment a new mild carbon flash is occuring rela- 
tively close to the carbon helium discontinuity. The region 
in which carbon is burnt is located in the most external 
regions of the partially degenerate core, at only 0.3 Mq 
below the carbon helium discontinuity. Consequently, the 
layers between the carbon-helium discontinuity (which is 
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Fig. 8. Temporal evolution of the first 10 thermal pulses. 
The upper panel displays the luminosity associated to he- 
lium burning, wlicTC'as lower panel represents the evolution 
of the luminosity associated to hydrogen burning. 

located at Mr ~ 1.1936 Mq) and the BCE start expanding 
and cooling. Hence, the radiative gradient in this region 
increases over the adiabatic gradient and, consequently, 
the BCE advances towards the center. However, and op- 
posite to what occured in the previous carbon flashes, the 
third flash is not strong enough to halt the expansion and 
cooling of the region between the HeBS and the BCE, and 
therefore, the convective zone keeps advancing. 

As it has been already mentioned, after the second 
carbon flash, a third flash of smaller strength occurs. 
The maximum luminosity is in this case of about Lq « 
3.0 X 10^ Lq. Moreover, in this case the associated convec- 
tive zone develops in a region of the core which is closer to 
the BCE and, hence, mildly degenerate. After this flash, 
the remaining carbon burns out more quietly. The sub- 
sequent small decrease in the surface radius of the star 
(see the central panel of Fig. 5) and the overall heating 
of its envelope, as well as the presence of CNO elements 
in regions just below the base of the convective envelope 
allows a reignition of the hydrogen burning shell. The con- 
ditions of degeneracy allow for a hydrogen shell flash — 
see Fig. 8. This flash has a magnitude of Lh — 10^ Lq, 
lasts for 60 yr, and occurs at t ~ 8.065176 x 10^^ s in a 
region (M^ ~ 1.192763 Mq) just below the BCE (which 
is located at Mr ~ 1.192843 Mq), and very close to the 
locations of the almost extinguished carbon burning shell 
at Mr c:^ 1.178429 Mq and of the helium burning shefl, 
at Mr ~ 1.191962 Mq. The associated convective region 
is very small (~ 10~^ Mq) and rapidly merges with the 
convective envelope. At this time the carbon luminos- 
ity is Lc — 1.3819 X 10^ Lq and the helium luminosity 
is Lhc — 1.1071 X 10'* Lq. This flash is a direct conse- 
quence of the the compression of the layers just below the 
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Fig. 9. Abundance proflles of the core at the end of the 
carbon burning phase. 



BCE which, in turn, is a consequence of the extinction 
of the carbon burning shell. However, after the resurrec- 
tion of the hydrogen burning shell degeneracy is quickly 
removed due to the expansion and cooling which derive 
from the fast injection of energy in the layers involved. 
Afterwards, hydrogen burning is not quenched, but in- 
stead it proceeds at approximately constant luminosity of 
about Lh ~ 3.0 x 10** Lq. The reactivation of the hy- 
drogen burning shell increases the mass of the He buff'er 
and, consequently, marks the beginning of the alternance 
of hydrogen and helium shell flashes, which interchange 
themselves as main nuclear energy suppliers of the star. 
Therefore, our model star begins the thermally pulsing 
Super- Asymptotic Giant Branch, which will be described 
in detail in §4. The term Super- AGB stars was coined by 
Ritossa ct al. (1995) to refer to AGB stars with ONe cores. 
To this regard it is worth having a look at the abundance 
profiles of the inner ONe core. The core abundances after 
the main carbon-burning phase are shown in Fig. 9. The 
abundance pattern does not present striking differences 
with respect to the solar metallicity corresponding core. 
However, it is worth mentioning that the core of the zero 
metallicity star is slightly overabundant in ^°Ne with re- 
spect to the solar metallicity one — X(^°Ne)= 0.28 versus 
X(^''Nc)= 0.24 and somewhat underabundant in oxy- 
gen. Finally, the final mass of the C-exhausted core in the 
Z = case is Mc = 1.15 Mq, 0.09Mq larger than that of 
its Z = 0.02 counterpart. 

4. The thermal pulses 

Fig. 8 shows the temporal evolution of the hydrogen lu- 
minosity. We have computed 14 thermal pulses of which 

only the first 10 thermal pulses are shown in Fig. 8. After 
the first hydrogen fiash, a quiet phase during which he- 
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lium accumulates follows. During this phase the luminos- 
ity of the hydrogen burning shell remains nearly constant 
(Lh — 3.0 X 10^ I/q) and the helium burning shell is 
still active, opposite to what was found by Ritossa et al. 
(1995) for a solar metallicity star of 10 Mq. However, he- 
lium is radiatively burnt at the base of the helium buffer 
at a rate smaller than that of hydrogen at the base of 
its burning shell, the net result being that the mass of 
the helium buffer increases. After some time, the helium 
buffer reaches the critical mass and, consequently, a se- 
ries of thermal pulses begins. The first few pulses are 
rather weak, as usually found for this kind of stars — 
Ritossa et al. (1995) and Garcia-Berro et al. (1997) found 
some mini-pulses, which are totally equivalent to what 
is found here — but, as time advances, the amplitude 
of the pulses increases. Note as well that the inter-pulse 
time interval also varies from one cycle to the next. By 
the end of the calculations shown hcTc the amplitude of 
the pulses remained almost constant with an approximate 
value of i/He^ — 6-0 x 10^ Lq for the maximum luminosity 
and of about ^hc" — ^'^'^ the minimum luminos- 

ity, whereas the inter-pulse time interval did not, ranging 
between r » 250 yr and 400 yr. For the sake of com- 
parison, it is interesting to recall that Garcfa-Berro ct al. 
(1997) found a max;imum luminosity of about 6.0 x 10'^ Lq 
for their solar metallicity 9 Mq model star, and a much 
more regular inter-pulse time interval of 400 yr, whereas 
Ritossa et al. (1995) found a maximum luminosity of 
about 2.7 X 10^ Lq and an inter-pulse period of roughly 
200 yr for their solar metallicity IOMq model. Thus the 
thermal pulses of our Z = model star of 9 Mq are simi- 
lar to those of a model star of slightly larger mass of solar 
metallicity, as we have also found for other characteristics 
(coupling of the interior and the envelope during the car- 
bon burning phase, size of the core, chemical composition 
of the He-exhausted core. . . ) of our model star. 

The evolution of the temperature at the base of the 
convective envelope during the pulses is also midway to 
those of 10 Mq and 9 Mq model stars of solar composi- 
tion. During the mini-pulses it adopts a typical value of 
Tbce — 3 X lO'' K monotonically increasing to Tbce — 
8 X 10^ K up to the cigth pulse. From this thermal pulse 
on, the temperature at the base of the convective enve- 
lope remains roughly constant and enough to transform 
part of the freshly synthesized CNO clenic;nts into ^''Mg 
and neutrons. As to the existence of a third dredge-up, we 
have not detected such a process after following 14 helium 
flashes. However, it is important to realize at this point of 
the discussion that dredge-up processes are strongly de- 
pendent on the treatment of convection, and the fact that 
wc have not introduced any overshooting in our calcula- 
tions might have prevented the development of an addi- 
tional dredge-up. 

It is worth recalling at this point of the discussion that 
the thermally pulsing AGB phase of low and intermedi- 
ate mass primordial stars has been thoroughly studied by 
Siess ct al. (2002). They foimd two different mechanisms 
allowing for the thermal pulses. For objects between 1 



and 5Mq, a convective instability develops at the He-H 
discontinuity causing the ingestion of carbon and, con- 
sequently, reactivating the CNO-cycle. Instead, for the 
7 Mq object it is a previous dredge-up phase which al- 
lows the enrichment in CNO products and the subsequent 
behaviour as a normal thermally pulsing AGB star. Our 
9 Mq model star behaves very much in the same way as 
the 7 Mq star described by Siess et al. (2002). Both our 
results and those of Siess et al. (2002) contradict the pi- 
oneering work of Fujimoto et al. (1984). These authors 
found that there was a threshold value at = 0.73 Mq 
such that stars whose core mass was lower than M* experi- 
enced thermal pulses, but stars with a core more massive 
than this threshold value could only experience thermal 
pulses if their metallicity was above Z « 10^^. However, 
as already shown by Siess et al. (2002) and here, the pollu- 
tion in CNO-elements caused by convection in the He-H 
discontinuity zone or caused by previous dredge-up phases 
makes this requirement on initial metallicity unnecessary. 

Finally, it is interesting to have a look at the sur- 
face ratios of carbon, nitrogen and oxygen. For the solar 
metallicity stars of masses 9 and IOMq studied before 
in Ritossa et al. (1995) and Garcfa-Berro et al. (1997) 
these ratios are respectively (C:N:0)=(1. 63:3. 63:6. 93) and 
(2.35:4.25:6.26) after the second dredge-up, and the total 
metallicities are in both cases Z ~ 0.012 at this moment. 
In sharp contrast the initially Z — star of 9 Mq studied 
here presents (C:N:O)=(2.03:0.02:0.03) for a total metal- 
licity of Z = 2.1 X 10~^. Thus, observationally our model 
star would look like a C-enhanced metal-poor star. 

5. Mass loss and evolution to white dwarf 

In Fig. 10 the latest stages of the evolution in the 
Hertzprung- Russell diagram of the zero metallicity 9 Mq 
star are presented. The most important evolutionary 
stages have been highlighted along the evolutionary track. 
In particular, the beginning of the carbon burning phase 
has been marked with a circle, whereas the diamond shows 
the end of the thermally pulsing Super- AGB (TP-SAGB) 
phase. After the eleventh thermal pidse, when the sur- 
face radius is Rg ~ 320 -Rq and the surface luminosity is 
Lg ~ 6.5 X 10'^ Lq, we have assumed, somehow arbitrarily, 
that a mass-loss episode ensues. Before the eleventh pulse 
no mass-loss was assumed to occur. This mass-loss episode 
is clearly marked in Fig. 10 using a thick solid line. 

As it occurs for solar metallicity objects of analogous 
ZAMS mass, the mass-loss episode is plagued with many 
uncertainties, as no definite mass- loss rates are available 
for low metallicity stars, and not even there is a well stud- 
ied mechanism responsible for these mass-loss episodes — 
see, however, Marshall et al. (2004). Stellar winds for the 
case of single stars as well as Roche lobe overflow for the 
case of binaries are likely to play important roles, as it usu- 
ally occurs for solar metallicity stars of analogous masses. 
With the values of the luminosity and the radius at the 
end of the eleventh thermal pulse the mass-loss rates due 
to stellar winds would be significant for a solar metallicity 
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the early horizontal track corresponds to the mass loss 
episode. See text for further details. 



star, of the order of 10^*' Mq yr^^, according to Rcimcrs 
(1975) and Nieuwenhuijzen & de Jager (1990). We recall 
that the convective envelope of our model star presents a 
significant enrichment in metals {Z = 2.1 x 10~^) which, 
nevertheless, is smaller than that of a solar metallicity star 
of the same mass. Consequently, to a first approximation 
wc can consider the mass-loss rates of Nieuwenhuijzen & 
de Jager (1990) as representatives of the real mass-loss 
rate. Additionally, the extension of the mass-loss episode 
must also be taken with some caution and, consequently, 
the final limit of the thick line in Fig. 10 must also be 
considered as an approximate value. 

Keeping this in mind, we have performed different 
sets of calculations of the latest evolutionary stages, to 
cover the effects of mass loss with rates between ~ 
10-'' Mq yr-i and ~ IQ-* Mq yr^^ which bracket the 
representative value derived above. The evolutionary se- 
quences computed with low rates would describe the ef- 
fects of mass loss due to stellar winds, and the sets of evo- 
lutionary sequences computed with higher rates (between 
10~^ Mq yr~-^ and ~ 10~^ Mq yr~^) might describe mass 
loss due to Roche lobe overflow in a binary system. The 
mass-loss rate used for the evolutionary sequence shown 
in Fig. 10 is M = 3.0 x lO'^ Mq yr'^ Finally, it is im- 
portant to say that, in all the cases studied here we have 
encountered some computational difficulties at the latest 
evolutionary stages — as it is usually found for AGE stars 
with carbon-oxygen cores with masses smaller than those 
considered here — but no significant differences regarding 
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licity 9 Mq remnant star. 



the final size and composition of the remnant core have 
been found. 

Finally, Fig. 11 shows the chemical profiles of the rem- 
nant of the zero metallicity 9 Mq star. The mass of the 
core is ~ 1.19 Mq. If we compare with the central re- 
gions of the model presented in Fig. 9, corresponding to 
a time shortly after the end of the carbon burning phase, 
it can be seen that some sharp features in the abundance 
profiles after the bulk of the carbon burning phase have 
been smeared out due to the Ray leigh- Taylor instabilities 
(Salaris et al., 1998) that occur in the subsequent evolu- 
tion of the star, but neither the size nor the general com- 
position of the core have substantially changed. The main 
isotopes present in the core are ^^O and ^°Ne, with mass 
abundances ^ 0.59 and ^ 0.28, respectively. As in the 
case of the solar metallicity stars of similar mass (Ritossa 
et al., 1995; Gil-Pons & Garcfa-Berro, 2001), some small 
amount of unburnt carbon is present (with a maximum 
mass abundance of ~ 0.02 before the Rayleigh- Taylor in- 
stabilities mix it through out the inner regions of the core). 
Finally, in the Z = case no significant amounts of ^^Mg 
nor ^^Na have been found. 



6. Summary and Conclusions 

We have computed the evolution of a zero metallicity 
9Mq star from the zero age main sequence until the 
formation of an ONe white dwarf, through core carbon 
burning and the TP-SAGB phase. Previous works on zero 
metallicity intermediate mass evolution have recently been 
made by Marigo et al. (2001), Marigo (2002), Heger & 
Woosley (2002) and Chieffi et al. (2001), but these au- 
thors did not follow carbon ignition in the degenerate core. 
Hence, this is the first calculation in which the carbon 
burning phase and the thermal pulses of a heavy-weight 
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intermediate mass primordial star has been followed self- 
consistently. 

Our results compare favorably with those of Marigo 
et al. (2001) and Marigo (2002) for a 9Mq model star. 
In particular, the time scales for hydrogen burning are 
notably similar in both calculations. We do not obtain 
such a nice agreement for the final size of the CO cores, 
since our calculations yield a ~ 1.1 Mq CO core, whereas 
they obtain a CO core of 1.28 Mq. However, it should 
be noted that even though wc both are using the OPAL 
radiative opacities (Rogers & Iglesias, 1992; Iglcsias & 
Rogers, 1993), and the conductive opacities of Hubbard & 
Lampe (1969), our respective equations of state are not ex- 
actly the same. Concerning the nuclear reaction rates, we 
both arc using basically the reaction rates of Caughlan & 
Fowler (1988), but Marigo ct al. (2001) increased the rate 
for the controversial ^^C(q!, 7)^^0 reaction by a factor 1.7 
(Weaver & Woosley, 1993), whereas we did not. Finally, 
the factor that is more likely to ultimately determine the 
differences in the sizes of our respective CO cores, is proba- 
bly the treatment of convection. Marigo et al. (2001) used 
moderate convective core overshooting, whereas we did 
not use it at all. 

We obtain a partially degenerate CO core in which 
carbon is ignited off-center. Carbon burning proceeds 
through a series of flashes which generates inner convec- 
tive regions, as it is found for stars of similar mass of so- 
lar metallicity. Regarding convection, it is worth noticing 
at this point that a calibration of the parameters of the 
Mixing Length Theory (and in particular of the amount 
of overshooting) can only be made for the case of solar 
metallicity by comparing with the abundant observational 
material. This procedure provides reliability to the MLT 
results, but the problem in the case oi Z = stars is 
the lack of observations and, therefore, of an appropriate 
calibration. Given the lack of observational material we 
have preferred to stick as close as possible to our previ- 
ous calculations (Gil-Pons & Garcia-Berro, 2001; 2002) 
in which no overshooting was adopted. Although the evo- 
lution during carbon burning for the Z — model star is 
very similar to that of solar metallicity, it is important to 
note one difference between the two models. Namely, that 
the envelope of the Z = star behaves as if it were more 
decoupled from the core than its solar metallicity counter- 
part. Hence, its structure is not significantly changed by 
the release of energy that occurs at the innermost regions 
of the star. Once carbon is depleted in most of the central 
regions of the He-exhausted core we find that the base of 
the convective envelope advances inwards and dredges up 
nuclear processed elements. Finally, the hydrogen burning 
shell is reactivated and a thermally pulsing Super-AGB 
(an AGB star with an ONe core) is formed. We have fol- 
lowed the first 11 pulses and we have found that their 
characteristics arc midway between to those computed for 
solar metallicity stars of 9 and 10 Mq. 

Somehow arbitrarilly we have chosen mass loss to be- 
gin at t = 8.066307 x 10"'^^ s, when the eleventh thermal 
pulse is already over and the surface radius is approxi- 



mately 320 i?Q. At this stage, the dredge-up process fol- 
lowing core helium exhaustion has already taken place. We 
have not delved into the causes of mass loss, that might 
correspond either to stellar winds or to Roche lobe over- 
flow (in binary systems). Instead, because no parametriza- 
tion of these processes has been made up to now, we have 
computed a series of sequences with mass loss rates be- 
tween 10~^ Mq yr~^ and ~ 10~* Mq yr~^. In any case, 
no significant differences have been found regarding the 
mass of the remnant and its composition. We obtain a ul- 
tramassive helium-rich 1.19 Mq white dwarf formed by a 
1.17 Mq oxygen neon core and a 0.02 Mq carbon-oxygen 
buffer. Such ultramassive white dwarfs have been detected 
by Finley et al. (1997) in a spectroscopic survey of DA 
white dwarfs hotter than ^ 25000 K, and by Dupuis et 
al. (2002) in the EUVE and ROSAT surveys. More re- 
cently, Liebert et al. (2004) have found a tail of massive 
white dwarfs in the mass distribution of white dwarfs in 
the Palomar Green Survey. Moreover, these authors have 
pointed out that of 28 white dwarfs all but nine are within 
±30° of the Galactic plane, thus suggesting that there is 
an association between these white dwarfs and that of 
B and O stars. Additionally, they noticed that the ul- 
tramassive component of the white dwarf mass distribu- 
tion has small cooling ages, suggesting massive progeni- 
tors. In summary, all these surveys have provided a hand- 
ful of ultramassive white dwarfs and it is interesting to 
realize that the derived mass distributions show a high 
mass shoulder with a peak centered at ~ 0.9 Mq that 
provides an excess of white dwarfs with masses between 
1.0 and 1.2 Mq. Since the theoretical mass distributions 
using a steady initial-to-final mass relation do not predict 
such behavior, Weidemann (2000) was the first to ascribe 
this behavior to an unsteady initial-to-final mass relation 
which separates white dwarfs with carbon-oxygen cores 
from those with oxygen-neon cores like the one obtained 
in the present work. Although the unsteady behavior has 
also been ascribed to binary evolution, the matter still re- 
mains as an open question, as recent Smoothed Particle 
Hydrodynamics simulations (Guerrero et al., 2004) do not 
totally discard the formation of an ultramassive white 
dwarf. Our simulations, thus, reinforce the idea that ultra- 
massive white dwarfs could be the end-product of the stan- 
dard evolution of single heavy-weight intermediate mass 
stars. 
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